-8- 



REMARKS: 

The present Amendment makes editorial changes in the title, 
specification and claims, and adds an Abstract, to conform the present PCT 
application to the requirements of United States patent practice. The 
5 cancellation of claims 1-16 in favor of the claims presented herein has been 
undertaken solely because the amount of strikethroughs and underlining that 
would have been necessary to confirm the original claims to the requirements 
of 35 U.S.C. §112, second paragraph would have been unduly burdensome 
and confusing. No change in the claim language has been made for 
10 distinguishing any claim over the teachings of the prior art of record. 
Accordingly, no change difference in any claim language between the 
language of the original claims and the language of the claims presented 
herein is considered by the Applicants as a surrender of any of the subject 
matter encompassed within the scope of original claims. 

15 Early consideration on the merits is respectfully requested. 
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D e script i on 

Method and arrangomont and computor program with program codo moans and 
comput e r program product for the analys i s of n e urona l activit ie s i n n e uronal areas 

SPECIFICATION 

5 TITLE 

"METHOD, ARRANGEMENT AND COMPUTER PROGRAM PRODUCT FOR 
ANALYZING NEURONAL ACTIVITIES" 

BACKGROUND OF T HE INVENTION 

Field of the Invention 

10 The present invention relates to aft analysis of neuronal activities in 

neuronal areas, for example of nerve structures in areas of the brain of a patient. 

Description of the Prior Art 

Knowledge about a mode of operation of a neuronal area and about an 
interaction of neuronal areas are fundamental to functional nucl e ar sp i n magnetic 
1 5 resonance tomography (fMRI), as described in A. W. Toga and J. C. Maziotta (Ed.), 
"Brain Mapping: The Methods", Ch 9: M. S. Cohen: "Rapid MRI and Functional 
Applications", Academic Press 1996 or fMRI t e chnology [3] , which is a further 
development of th e known magnetic resonance tomography. 

Pr e v i ous l y known magn e t i c Magnetic resonance tomography (a l so nuc l ear 
20 sp i n tomography, abbr e v i at e d to MR) is an imaging method which that produces 
sectional images of the human body without using damaging X-rays. 



Instead MR takes advantage of the behavior of body tissue in a strong 
magnetic field. Pathological changes in body tissue, fof example in the brain or 
spinal cord, can be detected by this modality m e ans . 



25 Functional disorders in body tissue, particularly in the brain of a patient, 

however, cannot , how e v e r, be detected by means of conventional magnetic 
resonance tomography. 
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This is performed by functional nuc le ar — sp\& magnetic resonance 
tomography or fMRI technology. 

Tho neurona l Neuronal activity in areas of the brain of a patient can be 
measured indirectly by means of the fMRI technique. The BOLD (Blood 
5 Oxygenation Level Dependent) signal, as it is called, is measured in individual 
areas of the brain, sa\4 this signal relating to neuronal activity in the respective 
areas. 

Dependencies, which stem among other things from structures in the brain, 
i.e. from neuronal links between nerve cells or nerve structures, exist between the 
10 neuronal activities in the areas. 

The outcome of the fMRI measurements shows the course of activity of 
individual areas over a certain period, for example during cognitive sequences as 
a result of specific perception processes or motor tasks. 

Functional disorders, in this case in the brain, are thus imp li c i tly inherently 
15 contained in the fMRI signals measured. 

Efficient methods for the analysis and evaluation of such fMRI 
measurements are thus desirable in order to be able to furnish evidence of 
possibly existing functional disorders in specific areas. 

Pr e vious l y known Known methods, such as described in A.R. Mcintosh et 
20 al., Structural Equation Modeling and Its Application to Network Analysis in 
Functional Brain Imaging, Human Brain Mapping, 2:2-22, 1994 for e xampl e th e 
m e thod of analysis known from [6] , are restricted to a detection of functional 
relationships between various areas of the brain in certain predetermined tasks 
such as the aforementioned perception processes or motor tasks (functional 
25 connectivity). These functional relationships are also designated functional 
connectivity. 
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In contrast to functional connectivity, however, the determination of a true 
physical connectivity, i.e. the determination of actually existing linking structures 
(of areas of the brain) independently of specific predetermined tasks, is not 
possible with these known methods. 

5 A further known method of analysis for detecting functional connectivity is 

described below. 

The object of this known method of analysis described below is the above- 
described detection of functional relationships between various areas of the brain 
in specific perception processes or motor tasks. 

10 This known method of analysis is based upon on a predefined model of a 

brain, i.e. a predefined brain architecture. 

This brain architecture, predetermined a priori from prior knowledge, 
defines general functional and/or spatial dependencies between specific areas of 
the brain in the form of a coupling matrix S, as it is called. 

15 The coupling matrix S has a form or structure that is fixed in accordance 

with the predetermined brain architecture (columns/rows) and is accordingly 
populated in certain, but not at all (matrix) positions^ with changeable coupling 
strengths Sj. These are changeable and are matched as part of the method of 
analysis. 

20 The unpopulated (matrix) positions are populated with fixed, unchangeable 

values, namely zero. 

The coupling strengths Si describe functional dependencies respectively 
between two areas of the brain or the BOLD signals measured there and 
r e pr e s e nt i ng represent the neuronal activities there. 

25 In this known method of analysis, the (changeable) coupling strengths Si 

are ftew defined such that statistical characteristic quantities wh i ch that are 
determined by this method of analysis from the fMRI measurements, can best be 
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explained. Expressed differently, a the probability for an occurrence of the 
measured data, i.e. the fMRI measurement or the BOLD signals, is to b e 
maximized by means of the desired coupling strengths Si. 

In this method of analysis a data point s=s t represents a totality of all BOLD 
5 signals s1, sN of the individual n areas at a point in time t or averaged over a 
time interval t(t=[1;T]). 

The fMRI measurement compris e s includes a large number of such data 
points for possibly differing perception processes and/or motor tasks, for which the 
corresponding BOLD signals were measured. 

10 In the this known method of analysis, it is «ew not the individual data points 

s1, s2, sJ x which are evaluated directly^ but statistical characteristic quantities 
which that emerge from these. 

For a statistical distribution of the data points s1, s2, sT it is assumed 
that it is described fully by a multivariant normal distribution, i.e. a statistical 
15 distribution of the first order, with a mean value m and a covariance £: 

P{s\/u i X) = — = ^ *e 2 (1) 

V2tt •|Z| 

For sufficiently long series of measurements, the occurrence of the 
individual data points si of s1, s2, sT can be viewed as statistically 
independent. 

20 The probability P=P(s1, sT|p, Z) for an occurrence of all measured data 

points s1 , sT can accordingly be written as: 

T 

P(s x ,...,s T \ ju,X)=YlP{s t \ JU,I,) = 

t=\ 

• e 2 (2) 
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Here, the unknown variables, the mean value \x and the covariance Z, 
depend exclusively on a(brain) model which that describes the measured data. 

The model assumes a linear statistical relationship between the individual BOLD 
5 signals: 

N 

*i = T, S ij s j + *i for * =1 ' N 
or 

s = Ss + s (3) 

where s d e scrib e s designates the external influence on the individual BOLD 
10 signals, like a sensory input of sensory cells on the examined areas of the brain. 

The influence variables si and sj on different examined areas i and j may 
can in this case be correlated throughout. 

The model parameters to be specified are consequently the coupling 
strengths Sj of the underlying coupling matrix S, the mean value jus of the external 
1 5 influence e and the covariance Ze of e. 

The mean value \x and the covariance £ depend on these: 

Z = 4S,ZJ (4) 

20 In the this known method of analysis the model parameters are then 

determined such that the probability P=P(s1, sT| given in (2) for the 
occurrence of the measured data is maximized. 

For this purpose, a known maximum likelihood estimation method 
(optimization) is applied such as described in T.W. Anderson, An Introduction to 
25 Multivariable Statistical Analysis, Chapter 3, John Wiley & Sons, Inc., New York, 
London, Sydney, 1994 . 
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Using the relationships (4) in (2), an expression that is dependent on the 
coupling strengths Si, the mean value us and the covariance Ze is obtained, which 
expression is maximized by the optimization. 

The optimization then leads to the desired coupling strengths Sj between 
5 the BOLD signals. 

These in turn then enable detection of functional relationships between 
various areas of the brain in specific perception processes or motor tasks 
(functional connectivity). 

The This known method of analysis d e scrib e d h e r e inabov e exhibits the 
10 disadvantage , how e v e r, that the measured fMRI signals are e xp li c a b le able to be 
modeled only insufficiently accurately or that the model is matchable only 
insufficiently accurately to the measured fMRI signals^ and consequently the mode 
of operation or interaction of neuronal areas is only insufficiently mappable. This 
shortcoming could possibly lead to incorrect conclusions being drawn made with 
15 regard to connective functionality. 

A software tool for an fMRI method of analysis, an "fmri.pro", is known from 
Specification for "fmri.pro" software relating to quantitative fMRI analysis, 
obtainable on 07.09.2001, under http://www.med. uni-uenchen.de/radin/html/ 
arbeitsqruppen/fmri/ccfmri.html {4}. A device for implementing the fMRI technique 
20 is known from Specification of fMRI - device, obtainable on 07.09.2001, under 
http://www.unipublic.unizh.ch/campus/uninews/2001/0147/fmri.html 

SUMMARY OF THE INVENTION 

The An object of the present invention is cons e qu e nt l y to i nd i cat e provide 
an improved method of analys i s for the analysis of neuronal activities. The 
25 improved method of analysis is intended to be able to better explain measured 
fMRI signals and th e r e by thus to better describe the mode of operation and 
interaction of neuronal areas than is the case with the abov e above-described 
known method of analysis. 
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This obj e ct i s achi e ved in th e m e thod and th e arrang e m e nt and i n th e 
computer program with program codo moans and i n tho comput e r program 
product for th e analysis of n e urona l act i v i ti e s i n n e urona l ar e as hav i ng th e 
f e atures according to th e r e sp e ctiv e i nd e p e nd e nt Claim. 

5 I n th e This object is achieved in accordance with the invention by a method 

for the analysis of neuronal activities in neuronal areas using signals describing 
the neuronal activities, the wherein signals are determined , wh e r e by with each 
signal d e scr i b e s describing the neuronal activity in one of the neuronal areas. 

A matchable coupling which that is described by using matchable coupling 
10 variables wh i ch that describe a statistical relationship between the matchably 
coupled signals forms the basis of all the signals, not rather than just some of 
them. 

Probabilities for an occurrence of the signals are determined, whereby a 
statistical distribution forms the basis of the occurrence of the signals. 

15 All matchable coupling variables are determined by an optimization of the 

probabilities and hence matched. 

The neuronal activities are analyzed using the matchable coupling 
variables. 

The arrangement in accordance with the invention for the analysis of 
20 neuronal activities in neuronal areas by the use of signals describing the neuronal 
activities has units stand i ng functiona ll y i n contact connected with one another 
a&4 which are configured such that 

- the signals can be determined, whereby each signal describes the neuronal 
activity in one of the neuronal areas, 
25 - a matchable coupling can form the basis of all, not just some of the signals, 
said matchable coupling being described by using matchable coupling 
variables which describe a statistical relationship between the matchably 
coupled signals, 
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- probabilities for an occurrence of the signals can be determined, whereby a 
statistical distribution forms the basis of the occurrence of the signals, 

- all matchable coupling variables can be determined by an optimization of the 
probabilities and hence matched, 

5 - the neuronal activities can be analyzed using matchable coupling variables. 

It is essential to the invention that the matchable coupling^ which is 
described by using the matchable coupling variables., forms the basis of all 
signals. By this means, absolutely all coupling variables are determined by the 
optimization of probabilities and hence matched. 

10 I n th i s way th e The invention thus differs from the known method of 

analysis described here i nabov e i n that above because in the known method ef 
ana l ys i s matchable statistical coupling forms the basis of only some of the signals. 

Only these can be determined by means of the optimization of probabilities 
and hence matched. 

15 Seen c l early Thus , the known method of analysis thus presupposes a 

known, predetermined and stipulated neuronal structure. 

In contrast to this, no predetermined and fixed coupling structures are 
assumed in advance by the inventive approach. These emerge only within the 
context of the optimization. 

20 Through the interaction of an optimization method and a search method, 

i.e. the search for existing couplings and the determination of their optimum 
values, in the inventive approach both the coupling structure most probable on the 
basis of the signals and a coupling strength of the specific couplings are 
determined. 



25 
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A particularly advantageous aspect of the inventive approach is that this 
approach is independent of other methods and of possibly defective prior 
knowledge. No prior knowledge, or just outline prior knowledge, of coupling 
structures is needed in the case of the invention in order to analyze the neuronal 
5 activities. 

As a result of the flexibility achievable by the invention in the matching of 
couplings, neuronal structures can be determined with greater precision and in 
greater detail. 

The computer program with program code m e ans according to the 
10 invention is equipped to perform all steps in accordance with the method of 
analysis accord i ng to of the invention if when the program is run on a computer. 

The computer program product with program code m e ans stored on a 
machine-readable medium is equipped to perform all steps in accordance with the 
method of analysis accord i ng to of the invention if when the program is run on a 
15 computer. 

The arrangement and the computer program with program code means 
which that are equipped to perform all the steps in accordance with the inventive 
method o f analys i s i f th e program i s run on a comput e r and th e comput e r program 
product with program code moans stored on a mach i ne-readab l e medium wh i ch is 
20 e quipp e d to p e rform a ll th e st e ps in accord a nc e w i th th e i nv e nt i v e m e thod of 
analysis if th e program is run on a comput e r are particularly suitable for 
performing the method of analysis according to the invention or one of its further 
developments explained h e r ei nb el ow below . 

Pr e f e rr e d furth e r d e v e lopm e nts of th e i nv e nt i on w il l e m e rg e from th e 
25 dopondont c l aims. 

The further d e v el opm e nts embodiments described h e r e inaft e r below relate 
both to the methods and to the arrangement. 
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The invention and the further d e v el opm e nts embodiments described 
h e r ei naft e r below can be implemented both in software and in hardware, for 
example using a special electrical circuit (ASIC) , 

Furthermore, implementation of the invention or of a including the further 
5 d e v el opm e nt embodiments described h e r e inaft e r is possible by means of a 
computer-readable storage medium on which is stored the computer program with 
program code m e ans which executes tho i nvent i on or furth e r d e v el opment . 

The invention , including the or oach further development embodiments 
described h e r ei naft e r below can also be implemented in a computer program 
10 product which has a storage medium on which is stored the computer program 
with program code m e ans which executes the invention or furth e r d e v e lopm e nt . 

The statistical distribution which that forms the basis of the occurrence of 
the signals can be a statistical distribution of a first or a higher order. The higher 
order can be achieved by using an Edgeworth expansion , as described in Samuel 
15 Kotz, Norman L. Johnson (Editors-ln-Chief), Cornish-Fisher and Edgeworth 
Expansions, Ch. 4, pages 188-192, Encyclopedia of Statistical Sciences, Volume 
2, John Wiley & Sons, 1982 {2} or a sum of normal distributions. 

In such a statistical distribution of a higher order not only do mean value 
and covariance have to be matched to comply with a quantity of data - as in the 
20 case of such a statistical distribution of the first order - but further higher-order 
parameters like moments and cumulants also have to be matched. 

It should be noted that the possibilities mentioned for achieving a higher 
order without restricting the generality are only two selected statistical 
distributions. Other possibilities are known to persons skilled in the art. 

25 In addition, in the sum of normal distributions, the individual normal 

distributions and thus indirectly the neuronal activities can be weighted. 
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Optimization can also be implemented by means of a maximum likelihood 
estimation , as described in T.W. Anderson, An Introduction to Multivariable 
Statistical Analysis, Chapter 3, John Wiley & Sons, Inc., New York, London, 
Sydney, 1994 {4} method. 

5 In the optimization, a relationship between the linear statistical relationship 

and the statistical distribution can be taken into consideration as an auxiliary 
condition. 

It is also expedient, because by this m e ans the biological model of 
neuronal structures can be mapped in a more real way, for external influences on 
10 the signals to be taken into consideration in the linear statistical relationship. 
Such external influences may, for example, be sensory inputs of sensory cells on 
the examined areas. 

The determination of signals in the invention, for example of BOLD signals, 
can be carried out by measuring signals or else by transferring and/or reading in 
15 existing signals. 

The invention and furth e r d e v el opment d e scr i b e d ar e is particularly suitable 
for using use in an fMRI technique, which is considerably improved and more 
powerful as a result. 

In the context of such an fMRI use or fMRI examination, the neuronal areas 
20 are areas of the brain with corresponding nerve structures of patients to be 
examined and diagnosed. 

In the fMRI examination using the inventive approach, BOLD signals are 
measured in various areas of the brain of a patient for defined perceptory or motor 
tasks carried out by the patient, which BOLD signals describe or represent the 
25 neuronal activities in the respective areas of the brain. These are evaluated or 
analyzed, whereby the signal coupling variables are determined. 
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Using the analysis results, in particular the signal coupling variables, 
functional as well as physical dependencies between areas of the brain can be 
detected and determined. These can further be used for a diagnosis concerning a 
functional disorder in an area of the brain of a patient, for example by comparing 
5 "disordered" dependencies with those of healthy persons. 

An e x e mplary e mbodiment of th e i nv e nt i on i s e xp l ain e d b e low and shown 
i n the F i gur e s, i n which 

F i gur e 1 shows a d e vic e for carrying out an fMRI accord i ng to an e x e mplary 

ombodimont, 

10 F i gur e 2 shows a sk e tch contain i ng st e ps i n a m e thod for ana l yzing BOLD 

s i gnals according to an exemp l ary ombodimont. 
Ex e mplary e mbodim e nt: Functiona l nucl e ar spin tomography (fMRI) 

DESCRIPTION OF THE DRAWINGS 

Figure 1 illustrates a functional magnetic resonance imaging apparatus for 
15 conducting an fMRI examination in accordance with the invention. 

Figure 2 is a flowchart for analyzing BOLD signals according to an 
exemplary embodiment of the invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The invention will be explained in the context of an embodiment for 
20 functional magnetic resonance imaging (fMRI). 

Fig. 1 shows a dev i ce an apparatus 100 for carrying out functional nucl e ar 
spin tomogr a phy or magnetic resonance tomography ( abbrev i ated to fMRI) and a 
functional nucl ea r sp i n tomograph or magn e tic r e sonanc e tomograph 100 . 

Basic principles of fMRI technology, which is a further development of the 
25 known modality of magnetic resonance tomography, are known from A. W. Toga 
and J. C. Maziotta (Ed.), "Brain Mapping: The Methods", Ch 9: M. S. Cohen: 
"Rapid MRI and Functional Applications", Academic Press 1996 {3}. 
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The nuc l oar spin tomograph apparatus 100 has a closed tunnel 110 which 
is incorporated in a magnet 1 20 m such a way that th i s magn e t it generates a 
strong magnetic field in the tunnel 110. 

The nucloar spin tomograph apparatus 100 also has a patient table 130 
5 that can be introduced into the tunnel 110, on which table a patient is placed for 
an examination. 

In addition, the nuc l ear spin tomograph apparatus 100 has a control unit 
131 wh i ch that enables monitoring and control of the patient table 130 in the 
examination, for example a controlled introduction of the patient table 130 into the 
10 tunnel 120. 

As further compon e nts, th e nuc le ar spin tomograph The tomography 
apparatus 100 also has a measuring device 140 for the measurement of BOLD 
(Blood Oxygenation Level Dependent) signals, an associated evaluation device 
141 for evaluating the measured BOLD signals, in this case a high-performance 
1 5 computer, as well as an operating or interaction interface device 142 for operating 
personnel and a display device 143 for displaying an examination result. 

The components of the nuc le ar spin tomograph apparatus 100 are 
functionally connected to one another, for example via signal lines or data lines 
150, via which the data and signals can be transferred. 

20 Based upon on the fMRI technique, the neuronal activity in areas of the 

brain of a patient can be measured and analyzed and a diagnosis can be made 
d e r i v e d from th i s by means of the funct i on al nucl e ar spin tomograph apparatus 
100 shown in Fig. 1 . 

To this end, the BOLD (Blood Oxygenation Level Dependent) signal in 
25 individual selected areas of the brain of the patient is measured by m e ans of the 
measuring device 140, &e\4 the BOLD signal being conn e ct e d with representation 
of the neuronal activity in the respective area. 
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The result of such fMRI measurements shows the curve of the activity of 
the individual areas over a certain period of time, for example during cognitive 
sequences as a result of specific perception processes or motor tasks which have 
to be carried out by the patient during an examination. 

5 Functional disorders in the brain of the patient are thus i mp l ic i tly inherently 

contained in the measured fMRI signals. 

Using the evaluation device 141, which makes available or implements a 
fvew the inventive method of analysis, the fMRI measurements, i.e. the BOLD 
signals measured in individual areas of the brain, are analyzed. 

10 This riew inventive method of analysis represents in th i s cas e an improv e d 

furth e r d e velopm e nt of an improvement over the known method of analysis 
described above. 

In the f*ew inventive method of analysis, brain activity is determined in the 
form of corresponding activation patterns in the examined areas in the brain 
15 and/or relationships between activation patterns in the examined areas and from 
that conclusions drawn directly of functional disorders in the brain and their 
causes. 

The ftew inventive of analysis made available by the evaluation device 140 
is based upon an extended and more flexible model of the brain, of the neuron 
20 structures in the brain and their behavior, in particular their interactions, on the 
basis of which the measured BOLD signal is analyzed and evaluated. 

The basic principles of the ftew inventive method of analysis and the model 
are explained below. 

The results or the conclusions of an examination are shown on the display 
25 device 143 and can be further processed by means of the operating and 
i nt e ract i on interface device 142 in connection with the evaluation device 141. 
They also serve as a basis for a medical diagnosis in respoct of a patient to be 
examined and diagnosed. 
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Basic principles of the ftew inventive method of analysis [[(]] as set forth in 
Fig. 2, steps 210 to 250, [[)]] 

It is po i nt e d out should be noted that the new inventive method of analysis 
is an improved further development of the e\4 known method of analysis 
5 described above. Consequently, i t is th e cas e b e low that - unless stated 
otherwise - the o\4 known and new inventive method of analysis match for thos e 
correspond in parts. Where matching parts are mentioned explicitly, they exhibit 
the labeling used previously hereinabove. 

Using the new inventive method of analysis 200 the fMRI measurements 
10 (210), i.e. the BOLD signals in examined areas of the brain of a patient, are 
analyzed (210 to 250) and/or compared with reference fMRI measurements. In 
this way, conclusions are drawn directly about functional disorders in the 
examined brain and their causes. 

The new inventive method of analysis 200, which generates statistical 
15 characteristic quantities, such as statistical correlations between fMRI 
measurements in various areas of the brain, is based upon on an extended and 
more flexible mathematical model (220) of the brain based upon the known 
mathematical model according to A. W. Toga and J. C. Maziotta (Ed.), "Brain 
Mapping: The Methods", Ch 9: M. S. Cohen: "Rapid MRI and Functional 
20 Applications", Academic Press 1996 (£4. 

In this extended model (220) of the f*ew inventive method of analysis, the 
coupling matrix S is populated in all (matrix)positions by changeable coupling 
strengths Sj. 

In the fvew inventive method of analysis 200 all - because they are also 
25 changeable - coupling strengths Si are determined in such a way that statistical 
characteristic quantities which are determined from the fMRI measurements can 
best be explained (210 to 250). 
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A data point s=s t represents the totality of all the BOLD signals s1, sN of 
the individual n examined areas at a point in time t (or averaged over a time 
interval t)(t=[1;T]). 

The fMRI measurement comprises a large number of such data points s1, 
5 s2, sT for differing perception processes and/or motor tasks for which the 
corresponding BOLD signals were measured. 

In contrast to the o\4 known method of analysis,, in which a multivariant 
standard distribution was assumed for the statistical distribution of the data points, 
in the new method of analysis 200 a weighted sum of normal distributions is 
10 assumed for the statistical distribution (220). 



equivalence of probabilities P=P(s|C1, CL, ja1 , jaL, £1, ZL ) (230) (cf. 
(2)) for the occurrence of measured data points s1, s2, sT depend on more or 
15 different parameters than the mean value jli and the covariance E of the eld 
known method of analysis. 

In the new inventive method of analysis 200 specific statistical variables, 
which can be calculated for the chosen statistical distribution, are now correlated 
with the model parameters, i.e. the coupling strengths Si, the mean value jis of 
20 external influence s and the covariance Ze of e. 

These include inter alia the means values ji1, jaL, the covariances X1, 
EL and all the moments and cumulants of the chosen distribution of a higher 
order. 



P(s|C 1 , J C Lj /i I ,..^ L ,Z 1 ,... s Z L ) = 



(5) 




In this case, the chosen statistical distribution and thus also the 
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An i mpl i c i t inherent relationship between the statistical distribution 
parameters and the model parameters to be determined emerges from this, in this 
case taking into consideration the distribution (5) and the extended model based 
upon the model according to (3). 

5 ju = {i(Ci,C L ,/u x ,...,ML,Z u ...,X L ) 
Z = Y,(Ci,....,C L ,/ui,...,/u L , J L l ,... J y Z L ) 

2 = Z(S,Z f Z) (6) 

10 In conformance with the e\4 known method of analysis, the optimum model 

parameters are now determined (240) in an analogous manner in the new method 
of analysis 200 using maximum likelihood estimation [1] through optimization or 
maximization of the probabilities (5). 

The basic principles of maximum likelihood estimation are described in 
15 T.W. Anderson, An Introduction to Multivariable Statistical Analysis, Chapter 3, 
John Wiley & Sons, Inc., New York, London, Sydney, 1994 {4}. 

The parameters to be taken into consideration for the optimization are the 
parameters of the chosen statistical distribution of a higher order, in this case the 
weighted sum of normal distributions, the desired model parameters and statistical 
20 variables, in this case the mean value m and the covariance Z from (6), via which 
the relationships between the model parameters and the statistical distribution (5) 
were produced. 

These relationships from (6) must be taken into consideration as auxiliary 
conditions in the optimization. 

25 The optimization then leads to the desired coupling strengths S s which that 

describe dependencies between the BOLD signals (250) and form the basis of the 
further evaluation and medical diagnosis (250). 
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An alternat i v e to th e e x e mplary e mbod i m e nt doscr i b e d i s i ndicated below. 

to As an alternative, in place of the weighted sum of normal distributions, 
the distribution of data points can aise be described by an Edgeworth expansion. 

The basic principles of the Edgeworth expansion are described in Samuel 
5 Kotz, Norman L. Johnson (Editors-ln-Chief), Cornish-Fisher and Edgeworth 
Expansions, Ch. 4, pages 188-192, Encyclopedia of Statistical Sciences, Volume 
2, John Wiley & Sons, 1982 fa. 

Although modifications and changes may be suggested by those skilled in 
the art, it is the invention of the inventors to embody within the patent warranted 
10 heron all changes and modifications as reasonably and properly come within the 
scope of their contribution to the art. 
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